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Abstract. In this work, we study the spectral properties of matrix Hamihonians gen- 
erated by hnearizing the nonhnear Schrodinger equation about soliton sohitions. By a 
numerically assisted proof, we show that there are no embedded eigenvalues for the three 
dimensional cubic equation. Though we focus on a proof of the 3d cubic problem, this 
work presents a new algorithm for verifying certain spectral properties needed to study 
soliton stability. 

Source code for verification of our comptuations, and for further experimentation, are 



available at http : //www . math . toronto . edu/ simpson/f iles/ spec_prop_code . tgz 
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1. Introduction 

The nonlinear Schrodinger equation (NLS) in M x M'^, 
(1.1) liJt + AV' + gm^)ij = 0, V^(0, x) = V^o(x), 

appears in many different contexts. In applications, it appears as a leading order approx- 
imation in nonlinear optics, many body quantum systems, and hydrodynamics. It is also 
intrinsically interesting as a canonical example of the competition between nonlinearity and 
dispersion. 

For appropriate choices of the nonlinearity : M — > M, the equation is known to possess 



soliton solutions, nonlinear bound states satisfying (1.1) with the ansatz 

^(t,x) = e^^*i?(x; A), 



where A > is the soliton parameter. It is conjectured that any solution of (1.1) with 
appropriate nonlinearity that does not disperse as t — )■ cxd must eventually converge to a 
finite sum of stable solitons. This is referred to as the "soliton resolution" conjecture, a 



notoriously difficult problem to formulate, see 36 . 

A natural property to investigate is the stability of the solitons. In 14,39,40 , a criterion 
for orbital stability is established. Briefly, it says that if the derivative of the norm 
of -R(x; A) taken with respect to A is positive, then the soliton is orbitally stable. The 
perturbation remains small in a particular norm, in the case of NLS. If this derivative 
is negative, the soliton is unstable. 

Though these results on the orbital stability of solitons are very powerful, relying on 
much of the variational structure of the equations, they have three weaknesses. The first 
is that they do not say if a perturbed soliton reaches an asymptotically constant state; 
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orbital stability only assures us that the perturbation remains small. The second is that this 
approach provides no information if the derivative of the norm vanishes. This is the case 
of the L?' critical focusing NLS equation, with g{s) = s'^l'^ and for saturated nonlinearities 
which possess minimal mass solitons. Finally, the orbital stability fundamentally depends 
on the underlying equation possessing a known variational structure. Though this is not a 



valid criticism for (1.1), it is a problem for other equations, such as those studied in 32 



Alternatively, results such as 3 4][8 16 24 25 , and many others, prove asymptotic stability 



of a soliton or a collection of solitons; the system converges to specific solitons to as t — )■ oo, 
and the rest of the mass disperses. Asymptotic stability is usually proven perturbatively. 
The leading order behavior of the perturbation to the soliton is governed by the linearized 
operator. First, linear stability is proven by assesing the spectrum of the linearized operator. 
Then the nonlinearity is shown to be dominated by the linear flow. The spectrum of the 



linearized operator of (1.1) with monomial nonlinearity was studied in [6]. 

Embedded eigenvalues of the linear operator are detrimental to proving the necessary 
linear estimates. Indeed, they obstruct the needed dispersive estimates, as demonstrated 
in [9]. Thus, it is standard to make the assumption that there are no eigenvalues embedded 
in the essential spectrum. It is known that such a condition cannot be proven directly 
using abstract properties of the linearized operator, but must in fact be directly related to 
algebraic properties of the soliton itself. In this work we develop an algorithm for studying 



the spectral properties of the operator appearing when one linearizes (1.1) about a soliton 
solution. Furthermore, we use this algorithm to prove the absence of embedded eigenvalues 
or resonances for four NLS problems. We collect these results in the following section. 

Remark 1.1. Though we only prove Theorem^for a small number of cases, our objective 
in this work is to present an approach for verfiying the spectral hypotheses required for 
soliton stability theory. 

1.1. Main Results. Our results hinge on a so called spectral property based on linearized 
matrix Schrodinger operators. The specific form, and its motivations, are developed Section 
[5j In general, this property can be formulated as: 

Definition 1.2 (The Generalized Spectral Property). Let d > 1. Given L± and a skew 
adjoint operator A, consider the two real Schrodinger operators 

C+ = -A + V+, £_ = -A + V_, 

defined by 

^+/ = ^[^+,A]/=^[L+A/-AL+/], 
C.f = ^-[L.,A]f =^-[L.Af - AL.f] 

and 

V± = ^x-Vl^±. 
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Let the real quadratic form for z = {u, vY £ x be 

13{z,z) = l3+{u,u) + B^{v,v) 
= {C+u,u) + {C_v,v) . 

The system is said to satisfy a spectral property on the suhspace lA C x if there 
exists a universal constant 5o > such that Vz G W, 



Biz,z)>6oJ (|Vz|2 + e-|yl|zp) dy. 



In this work, the skew adjoint operator is 

This has particular significance for the critical equation, though we employ it in super- 
critcal problems. This (mis) application is discussed in Section [7j 

Our results rely on the key observation of G. Perelman p^ll that 



JL 



Theorem 1. Given the JL operator, arising from the linearization of NLS about a soliton, 
where 

L_ 
-L, 



assume the L± operators satisfies the Spectral Property in the sense of Definition L2 Then 
JL has no embedded eigenvalues on the subspace U. 

Proof. Let us assume we have an embedded eigenstate Zem = (wcm, "^^em)^ ^ U corresponding 
to eigenvalue zTem, > Aq. Then, 

L—Vqjji iTem^emi 

Plugging directly into the form, 

^(^eiii) ^cm) (-^+^0111; ^cm) ~l~ (•^ — ^cni) ^em) 

= ^ {(A^cm,-^+^cm) + (-^+^em, AWem)} 

+ - {{AVem, L-Vem) + {L At;cm)} 

^ { ( AWenn ^em ) (^em; AWgui) -|- (fern; Altguj) ( AWenn ^em ) } 

= 0. 
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□ 

We remark that this holds not just for embedded eigenvalues, but for any purely imag- 
inary eigenvalue. Thus, if the spectral property holds, we are assured that there are no 
imaginary eigenvalues on the designated subspace. The subspace U will be set by our 
analysis of the spectrum in Section |4j 

Definition 1.3. Separately, we say that a linearized NLS problem satisfies of a spectral 
condition if it lacks both: 

• Embedded eigenvalues, 

• Endpoint resonances. 

Our second theorem, which relies on the first is: 



Theorem 2. The spectral condition holds for the 3d cubic equation, (1.1) with d = 3 and 
g{s) = s, linearized about the ground state soliton R. 



We adapt the methods of 13 to give a numerically assisted proof of this result. 

Remark 1.4. Though the main result of this paper will be to establish Theorem^ our al- 
gorithm can also be used to establish the spectral condition for the one dimensional equation 
with g{s) = s^'^ and g{s) = s^. Again, this is for the problem linearized about the ground 
state soliton. 

Separately, we establish that for 3d problems with nonlinearities satisfying the conditions 
necessary for the existence of a soliton, as discussed in [2], one need only test for embedded 
eigenvalues that are: 

• Near the endpoints of the essential spectrum, 

• On a sufficiently low spherical harmonic. 

In Appendix |X] we give a proof of the following result using positive commutator arguments 
otherwise known as Mourre estimates: 

Theorem 3. Given a Hamiltonian H, there exists some M > such that for \fi\ > M , 
there are no solutions such that 

Similarly, if d > 2, for any embedded eigenvalue, u^, there exists a K > such that the 
spherical harmonic decomposition 



^afc(r, ij,)(f)k{(j),e) 
consists only of harmonics with k < K. 

Remark 1.5. Such results are well-known using resolvent estimate techniques; however, 
our approach provides easily computable limits on M and K in terms of the soliton solution. 
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1.2. Organization of Results. In Sections [2| [3] and |4| we review fundamental properties 
of (1.1) and the associated linearized operator. 

In Sections 4.1.2 and 4.1.3, we collect results and adapt the techniques of 12 to prove 



properties of the discrete spectrum and the absence of embedded resonances. 



we prove an appropriate spectral property as in Definition 1.2 



G. Perelman's observation then rules out embedded eigenvalues. 



Finally, in Section 
based on the work in 
This proves Theorem 

In Appendices A.l| and A. 2 we use Mourre multipliers to eliminate large embedded 
eigenvalues and large spherical harmonics from the expansion of an embedded eigenvalue. 
Though these results have been known via resolvent methods for some time, we aim to 
collect as much analytic information about the spectrum as possible, providing bounds for 
future estimates and computations. An overview of our numerical methods with bench- 
marks is then presented in Appendix [Cj 

Acknowledgments. This paper is an extension of a result of a thesis done by the first 
author under the direction of Daniel Tataru at the University of California, Berkeley that 
arose from a discussion with Wilhelm Schlag and Galina Perelman at the Mathematisches 
Forschungsinstitut Oberwolfach. The first author is supported by an NSF Postdoctoral 
Fellowship. The second author is supported in part by NSERC. In addition, the authors 
wish to thank Gadi Fibich, Michael Weinstein, Ian Zwiers and especially Wilhelm Schlag 
for many helpful conversations throughout the development of the paper. 



2. Properties of the Nonlinear Schrodinger Equation 



In this section we briefiy review some important properties of (1.1). For additional 
details, we refer the reader to the texts [5 ,35 . 



In general, for nonlinearity (7 : M — )■ M, (1.1) possesses the following invariants for data 
ipo e and \x.\iIjo G L^: 



Conservation of Mass (or Charge): 
Conservation of Energy: 



where 



Git) 



g{s)ds. 
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Pseudo-Conformal Conservation Law: 

||(x + 2ztV)^||i2 -4^2 /" G{\u\^)d^=\\x^\\l2- [ e{s)ds, 

JR<* Jo 

where 

e{s)= [ {4{d + 2)Gm')-4dg{\^ljm\')d^. 



a 



Note that (x + 22)f:V) is the Hamilton flow of the linear Schrodinger equation, so the 
above identity relates how the solution to the nonhnear equation is effected by the 
linear flow. 

Detailed proofs of these conservation laws can be arrived at easily using energy estimates 
or Noether's Theorem, which relates conservation laws to symmetries of an equation. 

In this work, we restrict our attention to focusing nonlinearities, such that g{s) > for 
all s G M. These are the nonlinearities that can yield soliton solutions. Often, g{s) = s 
for some a > 0. We examine one instance of the power nonlinearity, the three dimensional 
cubic problem (a = 1). 

As noted, a soliton solution takes the form 

^(t,x) = e^^*i?(x; A), 

where A > and -R(x; A) is a positive, radially symmetric, exponentially decaying solution 
of the equation: 

(2.1) AR- XR + g{\R\'^)R = 0. 

For power nonlinearities, the existence and uniqueness of the ground state soliton is well 
known. Additionally, the scaling properties of this case permit us to take A = 1. We do 
this in all that follows. 

Existence of the soliton is proved by in p] by minimizing the functional 



T(V') = J iVT/'I'dx 
with respect to the constraint of flxed 

vw = y"[G(i^p)-^i^r]rfx. 

Then, using the minimizing sequence and Schwarz symmetrization, one obtains the exis- 
tence of the nonnegative, spherically symmetric, decreasing soliton solution. Uniqueness is 



established in 19 by ODE methods. 

An important relation is that Q{\) = Q{R{-; A)) and E{X) = E{R{-; A)) are differentiable 
with respect to A. This fact can be determined from the early works of Shatah, namely 



29 , 30 . By differentiating Equation (2.1), Q and E with respect to A, we have 



dxE = -XdxQ. 
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Variational techniques developed in [14] and [31] tell us that when (5(A) = E{X) + XQ{X) is 
convex, or S"{X) > 0, the soliton is orbitally stable. For S"{X) < the soliton is unstable to 
small perturbations. This stability (instability) directly is closely related to the eigenvalues 
of the matrix Hamiltonians resulting from linearizing NLS about a soliton. For a brief 
reference on this subject, see [34], Chapter 4. 



3. Linearization about a Soliton 

Let us write down the form of NLS linearized about a soliton solution. First, we assume 
we have a solution ip = e*'^*(-R + 0(x, t)). Inserting this into the equation, we have 

(3.1) ^(0)t + A(0) = -(?(i?2)0 - 2g'{R')R'Rei(l>) + 0(0^), 

by splitting up into its real and imaginary parts, then doing a Taylor Expansion. Hence, 
if (f) = u + iv, we get 

(3.2) 9, ^ JL [l 

where 

(-) -K-l'oX-°iJ)(V.° 

and 

(3.4) L_ = -A + X-V-, V.=g{R), 

(3.5) L+ = -A + A - y+, V+ = g{R) + 2g'{R^)R^. 

Alternatively, if we formulate the problem in terms of and 0*, 

(3.6) d, = iH 
where 

and 

(3.8) V, = g{R') + g'{R')R\ = g'{R')R\ 

The potentials in the two formulations are related by = Vi + V2 and = Vi — V2. 

There are many things we can immediately say about L_, L+, JL and "H. For a reference 
on the spectral theory involved, see Hislop-Sigal \15j or Reed-Simon [23j. First of all, both 
L_ and L+ are self-adjoint operators. Also, L_ is a non-negative definite operator and 
its null space is span{i?}. Note also that the functions ^ for j = 1,2, ...,d are in the 
null space of L_^.. By comparison with the operator A + A and using the fact that R 
decays exponentially, we see that the essential spectrum of Ti is the set (—00, A] U [A, 00) 
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from Weyl's Theorem, see [12] and [23]. Equivalently, the essential spectrum of JL is 
(— ioOjiA] U [i\,ioo). Indeed, a{JL) = iaiTi). Finally, using the fact that L_ is non- 
negative definite and looking at eigenvalues "H^, we see 

(3.9) L_L+u = 

However, this can be rewritten as 

1 1 
Tv = L'^L^Ltv = z/^f 

1 

for V = Ltu. Since the operator T is self-adjoint, we must have z/ G M U iM. 

Typically, asymptotic stability is studied with the following assumptions made on the 
matrix Hamiltonian: 

Definition 3.1. A Hamiltonian, is called admissible if the following hold: 

(1) There are no embedded eigenvalues in the essential spectrum. 

(2) The only real eigenvalue in [—A, A] is 0. 

(3) The values ±A are not resonances. 

Definition 3.2. Let NLS have nonlinearity g. We call g admissible at A if there exists 
a soliton, Rx, for NLS and the Hamiltonian, Ti, resulting from linearization about Rx is 



admissible in terms of Definition 3.1 



Remark 3.3. For simplicity in exposition, a matrix Hamiltonian, H, is said to be admis- 
sible if it satisfies several spectral conditions. One of these properties is a lack of embedded 
eigenvalues in the essential spectrum, which to establish we employ a spectral property as 
in 



131 



The spectral conditions of 1-i from Definition 3.1 are generally required to prove dispersive 



estimates for the evolution operator associated with the linearized Hamiltonian equation. 



which in turn are required to prove asymptotic stability results for solitons. See 24 , 25 
for further discussion. Let Pd and Pc be the projections onto the discrete and continuous 
spectrum of %. 

See Figure [T] for a description of the spectral decomposition for T-i resulting from lineariz- 
ing about solitons with subcritical, critical and supercritical stability properties. 

4. Spectral Properties of the Linearized Hamiltonian 

We now give more detailed and formal statements on the spectral properties of the 
operator under investigation. 



4.1. A Survey of Results on the Spectrum of l-L. 
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Orbitally Unstable Spectral Decomposition - Exponential Instability (Supercritical Behavior) 

Essential Spectrum 




Discrete Spectrum 



Figure 1. Plots of the spectral decomposition for l-L when Rx is exponen- 
tially unstable (supercritical behavior). 

4.1.1. An Analytic Result on the Spectrum ofH. We formalize the heuristic discussion from 
Section Isj with the following theorem from [12] . Let us write the operator as 



n = nQ + v 



A + A 









-V2 ' 





A + A 


+ 


V2 





In 12 , the authors proved the following properties of the spectrum: 



Theorem 4 (Erdogan-Schlag). Assume there are no embedded eigenvalues in the continu- 
ous spectrum of a{l-L). The essential spectrum ofH equals (— oo, — A] U [A,oo). Moreover, 
cr{l-L) = —aiTi) = o-{l-i) = a{l-L*) and cr{l-L) C M U zM. The discrete spectrum consists of 
eigenvalues {zj}^^^, < N < oo, of finite multiplicity. For each zj ^ 0, the algebraic and 
geometric multiplicities coincide and Ran(H — Zj) is closed. The zero eigenvalue has finite 
multiplicity. 

4.1.2. Absence of Embedded Resonances. This result is developed in the earlier work of 
Erdogan-Schlag and Agmon [I]. 

Define the space 
where 



i^'"^ = {/l|xr/eLn. 
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Then, we have the following Theorem, proved in 12 : 

Theorem 5 (Erdogan-Schlag). Let V\, Vi have sufficient decay at oo. Then for any /x such 
that \fi\ > X, [T-Lq — {fi ±iO))~^V : X_i_ — )■ X_i_ is a compact operator, and 

/ + (7/o-(/^±^0))-V 

is invertible on these spaces. 

The proof relies on a similar argument to a restriction theorem from harmonic analysis, 
which follows from a calculation using the specific structure of Ro{z) = (—A — z)~^. This 
strategy emulates closely that of the bootstrapping argument of Agmon, for scalar 
operators. 

4.1.3. Discrete Spectrum. We wish to show that the spectrum of H when linearized about 
the ground state has the discrete spectral decomposition Figure [T] 

Remark 4.1. For the 3d cubic nonlinearity, the structure of the discrete spectrum away 
from the essential spectrum has been verified numerically in [10], whose methods we recall 
briefly here. 



In 25 , using arguments derived from [22], it is shown that the discrete spectrum for 
supercritical exponents is determined by the discrete spectrum of L±. We present here 
a slightly stronger version that works for linearizations about a minimal mass soliton, 
R = -Rmin, in saturated nonlinearities. Though we will not numerically analyze any sat- 
urated nonlinearities in the current work, we present the generalized argument saturated 
nonlinearities are also of interest. 

The verification of the discrete spectrum heavily relies on the following result: following 
Theorem 6 (Schlag). Assume that L_ has no discrete eigenvalues on the interval (0, A] 



andT-i is a Hamiltonian as in (3.7) resulting from linearizing about a minimal mass soliton. 



Then, the only discrete eigenvalue for T-L in the interval [—A, A] is 0. 
Remark 4.2. A very similar theorem appeared in l25j proving the same result for the 



Hamiltonian H of the form (3.7) formed from g{s) = s with A = 1 and x G M^. The proof 
follows with minimal changes and is adaptable to many cases, hence we include it below for 
completeness. 

Proof. We argue by contradiction. To this end, assume H has an eigenvalue away from 0, 
say at E. Let A = 1 for simplicity. Then Ti^ has an eigenvalue at some value E"^ G (0, 1]. 
Hence, we have 

L^L^ue = E'^ue, 

for E"^ < 1. Since L_ is self-adjoint, we see that ue -L <f). By elliptic regularity, we have 
that Ue & -f^itc- Let P be the projection orthogonal to 0. Let A = PL^P. Using that 

ker(L+) = span{(9j0|l < j < d} 
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and 

(4.1) {n-z)-^ = {no-z)-'[i-Ui[i-U2J{no-z)-'Ui]-' 

xU2J{no-z)-'Ui]-'], 

we have 

ker(L+) = span{dxR, djR\l < j < d} 

since 

{R, dxR) = 0. 

Take Eq to be the unique negative eigenvalue for L+. Then, define 

g{a) = {{L+-a)-'R,R), 

which is well-defined and differentiable on {Eq, 1) since is orthogonal to the kernel of L+. 
We have 

g'{a) = - a)-^ R, R) > 

and 

g{0) = ^{R,dxR) = 0. 

Hence, g{0) = is the only for g in the interval {Eq, 1) since 

lim g{a) — ?■ — oo. 

a— S>i?o 

Conversely, if Af = af for some — oo < a < 1, a ^ and / G L^, then f -L R and 

{PL^P-X)f = {A-X)f = 0. 

Since 

EQ{f,f)<{L^f,f) = X{f,f), 

we know that X > Eq. If X = Eq, then / is a ground state of and hence not orthogonal 
to R. However, g{X) = 0, hence A = 0. So, A has a collection of eigenvalues at 0. Define 

Q = span{i?, Rj, Rx, ue}. 

We would like to show that dim(^) = d + 3. Since (p is orthogonal to all the other functions, 
we need only show that the equation 

d 

(4.2) CiUe + C2Rx + Cj+2Rj = 



has only the trivial solution cj = for all j. By applying L+ to (4.2), we see 

ciL+ue + C2R = 0. 
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Taking the inner product with ue-, we conclude Ci = 0. This imphes that C2 = 0. As a 
result, Cj = for j = 2, . . . ,d + 2. Now, if we can show that 

(4.3) sup {AfJ)<l, 

ll/llL2./Gg 

then by the Courant minimax principle, there would be at least d + 3 eigenvalues less than 
1 for A. However, we have shown there are exactly d + 2 of them. Note that neither 
the minimal eigenfunction for L+ nor (p itself are eigenvalues of A due to orthogonality 



arguments. Hence, if we can prove (4.3), we have proved the result. 
By our assumption on the spectrum of L_, 

{PL-JPf, f) < (/, /) 
for / 7^ 0. Since < 1 by assumption, we can prove the stronger result that 

{AfJ)<E^{PL-JPfJ) 
for all f = auE + b(j) + c ■ VR + dR\. To this end, we have 
{Af, f) = {L+{auE), auE + c ■ V0 + c/^a) 

+ (L+c ■ auE + c ■ VR + dRx) 
= E^{LZ^{auE), auE + c ■ VR + dRx) 

+ E^dRx, Lz\auE)) 
< E^{Lz\auE), auE + c-VR + dRx) + E^{dRx, Lz\auE)) 

+ E^Lz\c ■ VR + dRx), {c-VR + dRx)) 
<E^{PLz'PfJ) 

since LZ^ue -L V R and L_ is positive definite on ^ \ i?. □ 



In order to test the discrete spectral assumptions, we briefly recall the work of [10], which 
requires some numerical computation. First, let review the Birman-Schwinger method, let 
H = — X = —A — V for V > 0. Since we are looking for small, positive eigenvalues of 
L_, so take Hf = —o?f for < a < A so we have L_/ = (A — o?)f . Set U = y/V and 
g = Uf, then 

g = U{~A + a^)-'Ug. 
In other words, g & L"^ is an eigenfunction for 

K{a) = U{-A + a^)-^Ug 
with eigenvalue 1 where K > 0, compact. Conversely, ii g & satisfies K{a)g = g, then 

f = U--'g = i-A + a')-'UgeL^ 
and Hf = —o?j. The eigenvalues of K[ol) are seen to be strictly increasing as a — )■ since 

K\a) = -2XU{-A + a^)-^U. 
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This implies that 

# {a : KeT{H - a^) ^ {Q}} = # > 1 : Ker(ir(0) - E) ^ {0}} 

counted with multiphcity. 

Finally, use the symmetric resolvent identity to see 

{H-z)-^ = {-A- z)-^ + {-A- z)-^U 

x[l- U{-A - z)-'U] U{-A - z)-\ 

Hence, the Laurent expansion about z = does not require negative powers for z iff 
I + U{—A — z)~^U is invertible at 2; = 0, i.e. 

ker(/ - U{-A)-^U) = {0} 

by the Fredholm alternative since V has exponential decay. Hence, if H has no resonance or 
eigenvalue at 0, then K{0) will not have an eigenvalue at 1. If we then count the eigenvalues 
aj in decreasing order for K{0), then H has exactly negative eigenvalues and neither an 
eigenvalue or resonance at iff ai > 02 > ■ ■ ■ > «7V > 1 and un+i < 1- Hence, we can 
study numerically study K for a soliton of the saturated nonlinear Schrodinger equation. 
To do so, we must accurately find a soliton, then use it as potential in the truncation 
and discretization scheme presented in 10 , where the gap condition is verified for the 



Hamiltonian resulting from linearization about the 3d cubic ground state soliton. 

4.2. Generalized Kernel. Let us review the generalized kernel of a Hamiltonian resulting 



from linearizing about a soliton. Following 39 , we see by direct calculation that the vectors 



' " 




Rj 


R 








for all j = 1, . . . ,d are contained in ker( JL). Now, as Q{R\) is differentiable with respect 
to A, we have by a simple calculation that L^dxR = —R and L_{x(j)) = — 2V-R. Hence, the 
vectors 








" id,R),, ' 


XjR 








in the generalized null space of order 2. Notice that so far we have constructed at 2d + 2 
dimensional null space. Since we know the null spaces of L_ and exactly, these are 
unique. For power nonlinearities, g{s) = s'^, 



(4.4) 



(dxR) 



Ao=l 



-i- 

2V 



i? + x- Vi?). 



We use this explicit form in our calculations. 
As a result, we have the following 
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Theorem 7. Let g be the supercritical monomial nonlinearity. There exists a 2d + 2 
dimensional null space for Ti, the matrix Hamiltonian resulting from linearization about 
the ground state soliton (X = 1), consisting of the span of the vectors 








Rj 









R 









XjR 









Proof. The generahzed nuU space of the adjoint can be found by reversing the location of 
the non-zero elements in the above vectors. 

Suppose that there exists a generalized eigenspace for the eigenvalue E 0. Then, there 
exists X 7^ and i' such that ("H — E)ip = x and ("H — E)x = 0. Then, note that 



{A + E)x = 2Ex, 
0. 



Hence, has a generalized eigenspace at E'^. As a result, we see that T = L^L^ has a 
generalized eigenspace at E"^. Let Tx = E'^x and (T — E'^)ip = ex, for some c 7^ 0. Hence, 



{LlL+Ll - E^)Ll^i 



cXi, 



0, 



where given = I — Pr, we have xi = PrX 7^ since Tx = E'^x and ipi = P^ip 7^ 

1 1 

since (T — E'^)tp = cx- However, this means that the self-adjoint operator LlL^Ll has a 
generalized eigenvalue, which is impossible by an orthogonality argument. 

Since we have assumed there are no eigenvalues at the endpoints of the continuous 
spectrum, there can be no accumulation and the number of discrete eigenvalues is finite. □ 



4.3. Natural Orthogonality Conditions. As noted in Definition L2 even if the spectral 
property holds, it only implies the absence of embedded eigenvalues on a subspac e of 
U G L"^ X L^. That we must limit ourselves to a subspace will become clear in Section 5.2 



where we demonstrate that the operators C± have negative eigenvalues. 

This subspace will be defined as the orthogonal complement to the span of a set of 
vectors. If this collection of vectors is not chosen properly, we may find that the spectral 
property holds though the operator still has embedded eigenvalues. Thus the constraints 
on the set of vectors whose orthogonal complement will define U are: 

(1) They must be orthogonal to any embedded eigenvalues, 

(2) Orthogonality with respect to them should induce positivity of C on U. 

A way of meeting both of these requirements is to use the discrete spectrum of the adjoint 
matrix Hamiltonian, Ti*. To that end, we rely on the following simple results. 



16 J.L. MARZUOLA AND G. SIMPSON 

Lemma 4.3. If{X,u) is an eigenvalue, eigenvector pair for J L and {(t,v) is an eigenvalue, 
eigenvector pair for {J L)* , then 

(A - a*) {u,v} = 0. 
Thus, if X — a* ^ 0, the states are orthogonal. 

Corollary 4.4. An eigenstate of JL associated with an imaginary (possibly embedded) 
eigenvalue, ir ^ 0, is orthogonal to kerg((JL)*). 

Corollary 4.5. Let [ir ^ 0,-0) and (A > 0,0) be eigenvalue, eigenvector pairs of JL. Then 

(^1,02) =0, 
(^2,0l)=O. 

Proof. By the Hamiltonian symmetry of the problem, —A (02, 0i)"^ and A, (—02, —<PiY' are 
eigenvalue pairs of the adjoint, {JL)* . Therefore, 

(V'1,02) - (V^2,0l) = 0, 
-(^l,02)-(^2,0l)=O. 

Adding and subtracting these equations gives the result. □ 

These trivial observations motivate using the known spectrum of the adjoint system in 
constructing the orthogonal subspace. For the 3D cubic problem, we can thus use make 
use of eigenstates coming from the origin and the two off axis, real eigenvalues. 

5. Bilinear Forms and the Spectral Property 



We show here how the Spectral Property 1.2 is a condition sufficient for showing there 
are no embedded eigenvalues, thus proving Theorems [T] and [2j 

For general nonlinearities, we assume the operator resulting from linearizing about a 
soliton R has as discrete spectrum that is one of the following: 

• (i) a 2(i + 2 dimensional null space given by i?, V-R, d\R, x_R plus 2 eigenfunctions 
with symmetric discrete eigenvalues Aq, — Aq such that Aq G M, 

• {ii) a 2d + 4: dimensional null space given by R, VR, dxR, xi?, a, (3. 

In the following subsections, we establish the following: 

Theorem 8. The generalized spectral property holds for the 3d cubic problem for (/, G 
U <Z L^ X L^ specified by the following orthogonality conditions: 

(/,i?)=0, {g,R + ^-VR) = 0, (/,02)=O, {f,x,R) = forj = l,...d, 
where cj) = (0i,02)"^ is the eigenstate associated with the positive eigenvalue a > 0. 

This subspace is motivated by the observations in Section [4l3 as all of the elements we 
are orthogonal to arise from the spectrum of the adjoint problem. 
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Remark 5.1. Note that henceforward we assume here the unstable eigenfunction for the 3d 
cubic problem is radial. This claim is substantiated by direct integration in our numerical 
results section, as well as the fact that the spectral decomposition for the critical problem 
remains valid under the assumption of radial symmetry, hence the multiplicity of radial 
eigenfunctions must be at least 4. Since R, Rx are the only radial components of the 
kernel, the unstable eigenmode must also be radial. 



Following 13 , we computationally verify the spectral property in the following steps. 
First, the bilinear form is decomposed by spherical harmonics into 



(5-1) =J2Bi'\f^'\f^'^) + J2B['\9^'\9^% 

k=0 k=0 

where z = {f,g)^, and f^''^ and g^''^ are the components of / and g in the k-th spherical 
harmonic. We then identify the dimension of the subspace of negative eigenvalues for 
C^^\ Though at first this would appear to require an infinite number of computations, a 
monotonicity property of these operators with respect to k limits this to a finite number of 
harmonics. We then show that our orthogonality conditions are sufficient to point us away 
from the negative directions, allowing us to prove our result. 

5.1. The Index of an Operator. For a bilinear form 5 on a vector space V, the index 
of B with respect to V is given by 

indv'(i?) = max{A; G N | there exists a subspace P of codimension k 

such that -B|p is positive}. 

Our results rely on the following generalization of Theorem XIII. 8 of [23], which is in turn 
an extension of the Sturm Oscillation Theorem (Section XIII. 7 of [23]): 

Theorem 9. Let U^'^^ be the solution to 

dr^ r dr 

with initial conditions given by the limits 

lim ^ = 1, lim ^ = 0, 

r— >o r r— >o dr r 

where V is sufficiently smooth and decaying at oo. Then, the number N{U^''^) of zeros of 
U^^"^ is finite and 

indHi^j5(°)) = iV(?7(°)), 
ind^^i^^^(5('=)) = iV(t/(^)), A; > 1, 

where B'^^'> is the bilinear form associated to L^^\ 
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The space H^^^ is the set of radially symmetric H^{R'^) functions. The space -ffrad+ 



the subset of H^^^ for which 



X 



d'x < oo. 



We will omit the subscript notation in our subsequent index computations. It will be H^^^^ 
for = and -f^rad+ for > 1. 

If one wishes to remove the limits from the statement of the initial conditions, let 
U^^\r). Then the operator becomes 

_ _<P^ _ d-l + 2k d ^ ^ 
dr"^ r dr 

and the initial conditions become U^''\0) = 1 and ■^U^'^^O) = 0. Indeed, we use precisely 



this change of variables when making our numerical computations; see Appendix C.l The 
proof can be adapted from the proof of Theorem XIII. 8 in (23) . 

Corollary 5.2. The index is monotonia with respect to k, 

ind(5(^+i)) < ind(5('=)). 

This has the useful consequence that once we find an k for which ind(i?'''^^) = 0, we can 
immediately conclude that B^^ ^ > for all k' > k. 

Once we have computed the number of directions of each C^^"* that prevents it from being 
positive, we can check that we have a sufficient number of orthogonal conditions to point 
us into the positive subspace. 

5.2. Numerical Estimates of the Index. To compute the indexes of the operators, we 
proceed as follows. 

In dimension three, we solve the initial value problems 

(5.2) £ff/(o) = 0, f/(°)(0) = l, ^U^'\0) = 0, 

dr 

(5.3) £(°)Z(0) = 0, ZW(0) = 1, ^Z(°)(0) = 

dr 

for radially symmetric functions f/*-°^ and Z^^\ For higher harmonics. A; > 0, we solve the 
initial value problems 

(5.4) £f f/^'^) = 0, t/(^')(0) = 0, \imr-''U^''\r) = 1, 

r— >-0 

(5.5) = 0, Z^^\0) = 0, limr-'=Z('=)(r) = 1 
for radially symmetric functions f/*-'^-' and Z^''\ 
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Proposition 5.3. The indexes of 3d Cubic NLS are: 



ind/:L°^ = 1, ind£L^^ = 0. 



Once this proposition is established, Corollary 5.2 immediately gives us 
Corollary 5.4. For 3d Cubic NLS, 

mdB^_^^ =0 fork>2, 
ind-Si^^ =0 fork>l, 

where B^'' is the bilinear form associated with C^^^ . 

Using the method discussed in Section [c| we compute U^"^ and Z^"^ for each problem, 
a = 0, 1, 2. The profiles appear in Figures 2pp All were computed with a tolerance setting 

10-13. 

As a consistency check on the numerics, we note that asymptotically, the potential van- 
ishes, and 

(5.6) ^i^)^_^_d_^u^.)^Hk + d-2) 



In the region where r ^ 1, and the equations are essentially free and the solutions must 
behave as: 

(5.7a) U^'\r) ^ cl^\^ + cf V^-'^'-'^ 

and 

(5.7b) Z W (r) ^ r'^ + Z^j'^V^-'^-^ 

Estimating these constants from the numerics, we see that they have the "correct" signs. 
For instance in Figure [2] (a), f/^"-* clearly has one zero crossing. For sufficiently large r, 
the function appears to be increasing past a local minimum. However, since the constants 
appear to have stabilized, we contend we have entered the free region; the signs and mag- 
nitudes of the constants thus forbid another zero. 

Proposition 5.5. For the operators in Proposition |5.3| and Corollary 5^ there exists a 
universal 5q > 0, sufficiently small, such that for the perturbed operators 

the associated bilinear forms have the property that 

iYid{Bf) = ind(Sj^). 
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Figure 2. Index computations for 3d Cubic NLS. The number of zero cross- 
ings (other than r = 0), determines the codimension of the subspace on which 
the operator is positive. 




Figure 3. Index computations for 3d Cubic NLS. The number of zero cross- 
ings (other than r = 0), determines the codimension of the subspace on which 
the operator c'^^'' is positive. 



Proof. The proof here follows obviously from definition of the index of B, namely the 
positivity of the quadratic form B on the subspace for individual operators. Let 5q be a 

sufficiently small value that it holds for for k = 0, 1, 2 and for k = 0, 1. This 6q 
now holds for all higher values of k, again by a monotonicity argument. Indeed, for k > 2, 

E\'{fJ) = B\\fJ)+ / |/rcix>0. 
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Figure 4. Index computations for 3d Cubic NLS. The number of zero cross- 
ings (other than r = 0), determines the codimension of the subspace on which 
the operator d^^"* is positive. 



Because the form remains positive, this confirms that its index of zero is unperturbed. □ 



5.3. Invertibility of Operators. In conjunction with the results on the indexes of op- 
erators, we need to compute a number of inner products of the form {Cu,u), where C is 
one of our operators and u solves Cu = f. These are computed numerically, but we can 
rigorously justify the existence and unqiuess of these solutions, u, for the problems under 
consideration. 

Proposition 5.6 (Numerically Verified for 3d Problems). Let f be a smooth, radially sym- 
metric, localized function satisfying the bound \ f{r)\ < Ce"'^'" for some positive constants 
C and K. There exists a unique radially symmetric solution 

(1 + r'^+i)^ G L°"([0, cx))) n ^'([0, oo)) 

to 

(5.8) Cu = /, 

where C = for one of the 3d problems. 

Proof. This is Proposition 2 and 4 of |13j, along with our computations of the indexes in 
Proposition 5.3 See Appendix [B] for a proof in Id. □ 



Corollary 5.7. The solutions in Proposition 5j_6 are smooth and decay ocr as r ^ oo 
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5.4. Estimates of Inner Products. In order to prove the spectral property for each of 
these NLS equations, we need to approximate the bihnear forms associated with C^^^ on 
certain functions. These particular functions are, generically, of the form Cu = /, where / 
is from one of the orthogonality conditions. 

Proposition 5.8 (Numerical approximation of inner products). For the 3d cubic problem, 
let Uf, Uf \ U[^\ and zf^ solve 

(5.9) £f f/f ^ =R, (1 + r)f/f ^ G 

(5.10) £ff/f=02, (l + r)[/f gL°°, 

(5.11) di^U\''^ = rR, (1 + r2)f/f ) G L°°, 

(5.12) =R + rR', (1 + r)Z;°^ G 

where is the eigenstate associated with the positive real eigenvalue of JL. Then, 

(5.13) irf = (£f t/f ), f/f = 1.04846, 

(5.14) i^^f = (£f f/i°\ f/f ) = 0.00215981, 

(5.15) irf ^ (£f f/f f/f ) = -0.116369, 

(5.16) irf ^ = (^C^+^ui'\ f/f = -0.581854, 

(5.17) jf ^ = Zf^) = -0.662038. 

Proof. These result follows from direct computation. □ 
Finally, we state the following 



Proposition 5.9. For each case in Proposition 5^, there exists a Sq sufficiently small such 
that inner products associated with L can be made arbitrarily close to C These values will 
be denoted with overlines. 

Proof. This follows immediately from the invertibility of the operator and continuity. □ 

5.5. Proof of the Spectral Property. We are now ready to prove the spectral property. 
We prove positivity of the other cases are similar. Our proof closely follows Step 1 



and Step 3 of Section 2.4 of 13 . 

Since and K^"^ > 0, orthogonality to R and 02 will not give positivity. However, if 
/ is orthogonal to both of these, then it is also orthogonal to 
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and 



(0)^ A _ 7^(0) o^3 7^(0) , / ^3 r.{0) 



-5.22138. 



By Proposition |5.9[ we can take 5q sufficiently small such that 

1 ^1?^(0)n2 1>^(0)t.(0)^ 



^ 2 



((j?fy^-j?;°>j?f)<o. 



We proceed with this value of 5q. 
Let Q solve 

Obviously, 



and 



_ -(0) _ i^^^(O) 



For a moment, suppose Q G -ffradi not since it decays too slowly to be in L^. We 
could then imagine decomposing H^^^ into spanjQ} and its orthogonal complement, where 
the orthgonalization is done with respect to the Bj^ quadratic form. Since {Q,Q) < 0, 
the form is non- degenerate and this decomposition is well defined. Since indj^^i B^'f^ = 1, 

rad ~' 

-6+^ > on spanlQA}"*"- To prove this claim, we argue by contradiction. Suppose there were 

an element, Z G spanlQ^}"*" for which B^^\z, Z) < 0. Then, because of our decomposition, 

B^^^ < on span{Z, Q^}, a space of dimension 2. This contradicts our index calculation, 
proving the claim. 

Continuing, if m G -f^rad; -L (with respect to L^), then using the hypothetical orthog- 
onal decomposition, 

u = cQ + w^. 

If c = 0, then u lies in a subspace of H^^^ on which B^^^ > 0, giving the desired positivity. 
Indeed, the orthogonality condition, u -L q, is sufficient to ensure u is orthogonal to Q with 
respect to the B^^^ quadratic form. Taking the inner product of u with q, 



c (Q, q) + {u^, q) = cB';\Q, Q) + (n^, 

cBf (g, Q) + Bf {u^,Q) = cBf (g, g) + o. 
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Since Bf(Q,Q) ^ 0, we have c = 0. 

Unfortunately, the above argument does not work as stated because Q is not in To 

get positivity of B^^\ we regularize the problem and follow the above scheme. First, we 
introduce the smooth cutoff function XA{r) = x(r/A), defined such that 



1 r < 1 
r > 2 



and the norm 



II = liv/||i. + J |v± 

Let Q^(r) = Q{r)xA{r). Next, we observe that 

(5.18) hm \\Q^-Q\\+ + \S+\Q,Q)-'^+\Qa,Qa)\ = 0- 

Since Bf{Q,Q) < 0, for sufficiently large A, Bf{QA,QA) < too. Thus, we can 
legitimately decompose H^^^ as 

(5.19) i^rad = spanjQ^} © span{Q^}"^ 

with the orthgonalization is done with respect to the quadratic form B^^\ 

Finally, let u G H^^^. u ± R and ± 02- Then u A- q. With A sufficiently large to make 
the above decomposition valid, 

(5.20) u = c{A)Qj^ + ui. 
As in the heuristic argument b'^^\u\,u-^ > 0. Thus, 

Bfiu,u) = c{AfBf\Q^,Q^) + Bf{ui,ui) > c{AfBf{Q^,Q^). 

We will have our result if c{A) ^ as A — )■ +oo. 
Since {u, q) = 0, 

c{A) (U*^,q^^-{u\,q) 

I ± 7^(0) 7T 



Therefore, 
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Also, 

|(f7;,g) - {Q,q)\ = \(u*^-Q,cfQ)\ < iig -g^||+||g||+. 

Because this vanishes as A ^ +00, we have that for all A sufficiently large, 

\c{A)\<C\\ui\\4Q-Q4^ 

for a constant C independent of A. 
By construction, 

\\uiU<C + c{A)\\U*J+) < C + c{A)\\U*^ -QU+ WQU) ■ 

Substituting into our previous estimate on c{A), 

\c{A)\ < c {\\u\w + c{A)\\u\ -m+ + \\Q\U) \\Q - Qa\\+- 

We can clearly see that as A — )■ +00, c{A) — )■ +0. We conclude, 

Bf{u,u) > 

for u G H]^ and u -L R and m ± 02- This yields the estimate 
Following the same analysis for we conclude 

for g G i^rad g -L R + rR' since J^^^"* < 0. Repeating this again for / G -f^rad(i)' / 
we get 



because k']^^ < 0. 

Let us assume that 5q has been taken sufficiently small such that: 



The indexes of all operators are unperturbed, as in Proposition |5.5 
The above arguments on the positivity of for A; = 0, 1 and B^ hold. 
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Then for z = {f,g)'^ satisfying the orthogonahty conditions, 
B{z,z) = B4fJ) + B.{g,g) 



k=0 k=0 
oo oo „ 

k=0 k=0 

= Biz,z)-SoJ e-W'|z|2dx> 0. 

We almost have the expression in Definition L2 To complete the proof, note that for 
any (0,1), 

{l + 9)Biz,z)>9(^j\Vz\'d^ + Jv+\f\^ + j V^\gf^+5^ j e-\^\"\z\''di^. 
We can take 9 = 6^, sufficiently small such that 

0. + /v^-l^7r)+5oy"e-W>rdx>| |e-WVrdx. 



Then, 



Shrinking 5q again, so that it is smaller than 



mm 



gives us the spectral property. □ 

6. Other Problems 

In principle, this scheme can be applied to any linearized nonlinear Schrodinger equation. 
One finds the indexes of the operators, picks an appropriate subspace to project away 
from, and computes the necessary inner products. However, our experiments show that the 
algorithm is not as universal as might be hoped. In this section we exhibit the computations 
for several Id NLS equations, 

(6.1) i^t + i^.. + \^\'''i^ = o. 



Sometimes our approach works, ruling out embedded eigenvalues in a range of supercritical 
cases, while in others it fails, leaving a large range of interesting problems unresolved. 
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6.1. Numerical Estimates of the Index. As in the 3d problem, we first compute the 
indexes of the operators C± to identify the number of "bad" directions. In contrast to 
the multidimensional problems where there are an arbitrarily high, but finite, number of 
harmonics which must be examined, ID problems only require us to study the operators 
restricted to even and odd functions. This requires the following results, whose proofs are 
quite similar to that of Theorem [9] 

Corollary 6.1. Let U he the even solution to 

LU = -U" + V{r)U = 0, 
f/(0) = l,f/'(0) = 0, 

where V is sufficiently smooth and decaying at oo. Then, the number N{U) of zeros of U 
is finite and 

mdHiiB)=NiU), 
where B is the bilinear form associated with L. 

Corollary 6.2. Let U be the odd solution to 

LU = -U" - ^U' + V{r)U = 0, 



f/(0) = 0,f/'(0) = l, 

where V is sufficiently smooth and decaying at oo. Then, the number N{U) of zeros of U 
is finite and 

mdHi{B)=N{U), 
where B is the bilinear form associated with L. 

Hi and are the subspaces of if^(]R) restricted to even and odd functions. In what 
follows, we shall omit them in the subscripts of the indexes. 

To proceed, we numerically solve the initial value problems 

(6.2) = 0, U^'\Q) = 1, ^u^"\o) = 0, 

(6.3) = 0, Z^'\0) = 1, ^Z^'\0) = 

ax 

for even functions U^^^ and Z^^\ We then solve 

(6.4) £j^f/(°) = 0, U^"\0) = 0, ^U^°\0) = 1, 

ax 

(6.5) = 0, z^°\o) = 0, -^z(°)(o) = 1 

for odd functions U^°^ and Z^°\ C± have the same definitions as before; the Laplacian is 
now one dimensional. As in the 3d case, we verify a postiori that, asymptotically, [/" and 
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Figure 5. Index computations for critical Id NLS. The number of zero 
crossings (other than x = 0), determines the codimension of the subspace on 
which the operator >C^^°^ is positive. 



fit 
(6.6) 
(6.7) 

In addition, the constants have the right signs to ensure we have the correct number of 
zero crossings. 



Proposition 6.3 (Numerically Verified). The indexes for the Id NLS equation with a 
2,2.1,2.5,3 are: 



ind£j^ 
ind£i°^ 



ind£ 
ind£ 



0. 



Proof. Using the method discussed in Section [c| we compute U^"^ and Z^"^ for each prob- 
lem, a = 6,0. The profiles appear in Figures [7| [8j All are computed with a relative 
tolerance of 10~^° and an absolute tolerance of 10~^^ using Matlab . 

□ 



Proposition |5.5| applies to these Id problems too. As in the 3d case, we ultimately use a 
perturbed bilinear form in the proof of the spectral property. 



6.2. Estimates of the Inner Products. We now compute a series of inner products and 
show that in some cases the natural orthogonality conditions are sufficient to yield a spectral 
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Figure 6. Index computations for critical Id NLS. The number of zero 
crossings (other than x = 0), determines the codimension of the subspace on 
which the operator C^^^°'^ is positive. 
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Figure 7. Index computations for Id NLS with a = 2.1. The number of 
zero crossings (other than a; = 0), determines the codimension of the subspace 
on which the operator C^^^"^ is positive. 



property. Rigoursly, these results require the following Proposition on the invertibility of 
the C^^^°^ operators: 

Proposition 6.4 (Numerically Verified for Id Saturated Problem). Let f be a smooth, 
localized function satisfying the bound |/(a;)| < Ce"*^'^' for some positive constants C and 
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Figure 8. Index computations for Id NLS with a = 2.1. The number of 
zero crossings (other than a; = 0), determines the codimension of the subspace 
on which the operator C±.^"^ is positive. 




Figure 9. Index computations for Id NLS with a = 2.5. The number of 
zero crossings (other than x = 0), determines the codimension of the subspace 
on which the operator C^^^"^ is positive. 

K. If f is even/odd, there exists a unique even/odd solution u G L°°(M) fl C^(]R) to 
(6.8) Cu = /, 

where C = d'^^"'^ for a Id problem. 
Proof. See Appendix [Bj 
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Figure 10. Index computations for Id NLS with a = 2.5. The number of 
zero crossings (other than a; = 0), determines the codimension of the subspace 
on which the operator C^^^"^ is positive. 



a = 3 




(a) Even functions 
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Figure 11. Index computations for Id NLS with a = 3. The number of zero 
crossings (other than a; = 0), determines the codimension of the subspace on 
which the operator C^^^""^ is positive. 
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(a) Odd functions 
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(b) Odd function asymptotics 



Figure 12. Index computations for Id NLS with a = 3. The number of zero 
crossings (other than x = 0), determines the codimension of the subspace on 
which the operator C^^^"'^ is positive. 



Proposition 6.5 (Numerical). Let u['^\ z[^\ and Ui"\ all elements of L"^ 
following boundary value problems: 



solve 



(6.9) 

(6.10) 
(6.11) 
Let 
(6.12) 

(6.13) 
(6.14) 
Then, 



£(e)2-(e) 



R 

1 
a 



d 



dx 
R + xR' 

r(o) 



d 

dx 



f/r(o) = 0. 



£(e)^{e) 



a 



(o) 



2.0 -0.557768 0.292551 -1.30410 

2.1 -0.496932 0.216284 -1.21364 
2.5 -0.297841 -0.0216292 -0.924662 
3.0 -0.122559 -0.218499 -0.671783. 
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6.2.1. Proof of the Spectral Property for Certain Supercritical Cases. We restrict our at- 
tention to t he Id supercritical problems a = 2.5 and cr = 3. Repeating the procedure of 



Section 5.5, for z = {f,g)'^, the orthogonality of / to i? and xR gives us £+ > and the 
orthogonality of g to + xR! gives us £_ > 0. This proves the spectral property on the 
restricted subspace. Since these orthogonality conditions are consistent with those formu- 



lated in Section 4.3, we conclude that there are no non zero purely imaginary eigenvalues. 



6.2.2. An Inconclusive Supercritical Case. In the case of o" = 2.1, we have that jj'^ > 0, 
which means that orthogonality of g with respect to ^R + xR! , is insufficient to guarantee 
positivity of £_. It is possible that if we extend our scope, as in the 3D cubic problem, 
to include orthogonality to the eigenstate associated with the unstable eigenvalue we will 
be able to prove the spectral property for this problem. However, we do not pursue that 
here; rather we wish to highlight the failure of our algorithm at a seemingly arbitrary 
supercritical nonlinearity. 

6.2.3. The Critical Case. The critical ID problem, with a = 2, is also inconclusive. As in 
the supercritical problems we will look at the inner products against R and ^R + xR! . We 
also employ inner products arising from with the rest of the generalized kernel, x^R and /3, 
where /3 solves 

= -x^R. 



See 39 for details. This motivates the following numerical result: 



Proposition 


6.6 (Numerical). Let z'^'^ solve 


(6.15) 


&^zf =p, A^(^)(o) = 
ax 


and let 




(6.16) 


= B^^\Z^^\Z^^^) = i^&^z'^^\z'f 


(6.17) 


jf = B^1\Z^^\Z^^^) = l6^zf,zf 


Then, 




(6.18) 


= 3.77915, 


(6.19) 


4"^ = 0.864273. 



Since k[^^ < 0, we may conclude that £+ > 0, when the operator is restricted to even 
functions that are orthogonal to R. However, orthogonality to neither ^R + xR! nor p 

is, individually, sufficient to gain positivity of We are thus motivated to consider 
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orthogonality to the subspace span{|i? + xR' , p}, as in the proof of the 3D cubic problem. 
We examine the quantity 



(6.20) - ^ ((jf - ^ J^'^) = 0.0948958. 



J, 

However, we need this to be negative. Thus, we have no set of natural orthogonality 
conditions which yield a spectral property. In this case, since K'f'' < 0, the only obstacle 
to the proof is the C^f^ operator. 

6.2.4. The Critical Case with Other Orthogonality Conditions. If we had instead used the 
orthogonality condition, g R, and then solved the boundary value problem C^f^ z'^'^ = R, 
the inner product, 

j;^) = = -3.770731. 

This would give us a spectral property, but it is not a convenient subspace. 

Suppose we use the orthogonality conditions of 13 , and let g ± AR and g ± A^R. Then, 
we compute as follows: Let z['^\ and ^2^^ solve the following problems: 

(6.21) C^^^ZP = AR, -^Z['\0) = 0, 

dx 

(6.22) = AR, ^Zi'\0) = 0, 
and define the inner products 



(6.23) 




= B^^\Z['\ 




zi^^ 


(6.24) 


7(e) 
^2 


= B^^\zi'\ 


zf) = ^6^zf, 




(6.25) 


^3 


= B^^\Z^\ 


zf) = (^c^l^z^f. 




We will find that 










(6.26) 




fie) 


= 0.292551, 




(6.27) 




fie) 


= 2.57656, 




(6.28) 




fie) 
^3 


= -1.27657. 





As one would hope, given that the ID spectral property was established in 20 
(6.29) - ^ ((4'))2 - ) J^')) = -0.339932. 

This sign ensures that projection away from those two directions is sufficient to point us 

—(e) 

away from the negative eigenvalue, rendering £_ > 0. 
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7. Discussion 

We have demonstrated a computer assisted algorithm for proving the positivity of a 
bihnear form, B, on a subspace U. Because of the relationship between B, C, and the 
linearized operator, JL, we infer that there are no embedded eigenvalues. We succeeded 
with this program in the case of the 3d cubic equation, and a two supercritical Id problems. 
C. Sulem has suggested t that is likely to also be successful for solitons (with A = 1) of the 
3d cubic-quintic equation, 

for 7 sufficiently close to zero. We also anticipate success for other Id supercritical problems 
with a sufficiently large. These cases warrant further study. 

For subcritical problems, a similar algorithm should apply, though it will certainly require 
a additional orthogonality conditions. Many subcritical problems contian eigenvalues inside 
the spectral gap. Since our approach does not distinguish between embedded eigenvalues 
and imaginary eigenvalues in the gap, it would be essential to project away from those 
states. 

It remains to be seen how to extend our technique to other NLS/GP equations. Indeed, 
the failure in the ID critical problem is curious. The success or failure of the approach is 
likely related to the choice of our operator A = d/2 + x ■ V. In |T3], the authors proved 
the spectral property using this A, as it is generated by the scaling invariance of the mass 
critical problem. This results in the so-called "pseudoconformal invariant" for critical NLS 



and has great implications for blow-up. See 20 and 35 for additional details. 



Finally, recall that A determines the operators C±. These each have an index identi- 
fying the number of negative' directions. We then choose orthogonality conditions that 
simultaneously must satisfy the two properties: 

(1) They must be orthogonal to any embedded eigenvalues, 

(2) Orthogonality in with respect to these directions must imply orthogonality to 
the negative directions of C±, with respect to the B± quadratic form. 

The first requirement is satisfied by the vectors from the adjoint problem, as discussed in 



Section 4.3 We appear to have little fiexibility in altering these. Changing A will change 
C±; in turn this changes the negative directions. Thus, a different skew adjoint operator 
may extend the applicability of the algorithm. 
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Appendix A. Commutator Estimates 

A.l. Large Eigenvalues. We now establish upper bounds on the magnitude of embed- 
ded eigenvalues of JL. In this analysis, we examine the fourth order equation comes 
from squaring the operator. Therefore, we study eigenf unctions m G of the differential 
operator: 

(A.l) L^L+u = fx^u, 

where the L± operators are 

L± = -A + A - 14 

and /i G (Tcont.(<^-^) = (A,oo). As our proof applies to many NLS equations, we do not 
further specify the potentials V±; they are defined as in Section |3j 

From the properties of the soliton and the nonlinearity, we have that any solution to 
equation A.l is locally smooth via an iteration argument 18 . Following [2], asymptotic 



analysis shows that a solution decays exponentially fast. As a result, the possible range of 
frequencies is limited by 

(A. 2) ||Vu||i2 < yu||M||i2. 



See 33,37 for references on microlocal analysis. We prove the following: 



Theorem 10. There exists a fiQ > X such that for all fi > fiQ, the eigenvalue equation 
(A.l) has only the trivial zero solution in . 



Proof. We begin by defining the standard Mourre commutator as: 
(A.3) M = X ■ V + V ■ X. 

Using the structure of these operators, we immeadiately have the two identities 
(A.4) ([ML_L+ + L+L_M]u, u) = 0, 

(A.5) -/i>,M) = 0, 

where (■, ■) is the standard inner product. Combining these two identities with the 



frequency bound of (A.2), we can rule out solutions for sufficiently large fi. Indeed, 



((A' - 2AA + A(V+(x)) + V_(x)A- A(\/_ + V+) 



ML^L+ - L+L^M =[M,L_]L+ + L_[M,L+] + [L^,L+]M, 
[M,-A] = 4A, [M,\/_] = 2x- Vl^_. 

Hence, 

((-8A2 + 8AA - 4(\/_ + y+) A + 2(x ■ Vy_ + X ■ Vl^) A 

- 2A(x ■ VV- + X ■ VV+) + 2(x • VV^)V+ + 2(x ■ VV+)V- {**) 

- -y+)A- A(y_ -14)][d + 2x- V])m,m) = 0. 
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Since the last term is a product of skew-adjoint operators, they commute. 

Combining 4(*) + (**) and with the frequency bound ||Vit||i:2 < /i^||it||^2, we have for 
/i > Ho, this is a negative definite system. Thus, there are no eigenvalues. Furthermore, 
this estimate combined and standard Sobolev embeddings implies u e H'' for any k; hence, 
u is smooth. 

The system wc assess is: 



/ 



Hence, 



+ [4y_y+ - 2A(x • vy_ + x • w+) 

+ 2(x • vy_)y+ + 2(x • vy+)y_ + A(x • vy_ + x • vv+) 

- X • V(Ay_ - AV+) - d(AV^ - AV+)]u^ 
+ [4(x • Vii)(V(y_ - V+) ■ Vu) 

- 2(x • Vy_ + X • W+)S/u ■ Wu] dx = 0. 

-4{Auf - 4X{V^ + V+)u^)d^-4{fi^ - X'')\\u\\l2 + \\F\\Loo\\u\\l2 

+ (2||x- V\^_ +x- V^+IIloo 
+ Cama^\\dj{V^-V+)xk\\)\\Vu\\l, 



< J {-A{Auf - 4A(V_ + V+)u^)dyi 



-4(//2-A^-Ci-C2//')||ii"2 



L2) 



where 



F = 4T/_T/+-2A(x- vy_ + x- vy+) + 2(x- vy_)y+ 
+ 2(x • vy+)y_ + A(x • vy_ + X • w+) 
- X- v(Ay_ - Ay+) -d(Ay_ - Ay+) 

and Cj — Cj(y^, V^, A, d) for j — 1, 2. Hence, for ji large, we have that: 

j {-A{Auf - AX{V_ + V+)u'^ - Csiii, A, V_, V+)u^) rfx < 0, 

for C3 > and u a smooth function, hence u — 0. □ 

A. 2. Spherical Harmonics. As our potentials are radially symmetric, we expand our 
functions in spherical harmonics. Seperating the radial variable, r, from the angular vari- 
ables, 9, the expansion takes the form 
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where 



See 38 for a description of the eigenspaces of the spherical Laplacian, A5. Then, we have 



the following ODE eigenvalue problem: 



(A.6) 



92 
dr 

dr 



d-l d , / N 

d-l d 



X 



a 



r dr 



+ - + A-y+(r) -^^ 



Uk{r) = 0, 



where a"^ = k"^ + {d — l)d for k = 0,1,2, ... . Note, we use the notation since all of the 
eigenvalues are non-negative, which will be important in the sequel. 

We have following theorem: 



Theorem 11. There exists some ao > such that for all a > Oq, the eigenvalue equation 
(A.6) has only the trivial solution in L^. 



Proof. Let us denote the radial inner product by: 



and the operators: 



{u,v), 



uvr'^ ^dr. 



Ar 



dr 



dr 



d + 2r 



d_ 
dr 



Using the same commutator approach as in Section A.l 

4^2 



([4A,, 



-2rVL{r)][-/\r + 



a 



A- W r 



2 A 1 

+ [-A, + ^ + A-Vl(r)][4A,-^-2rV|(r)] 

+ [{V_ - F+) A, - A,(y_ - [d + 2r^]uk, Uk)r = 0. 
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Thus, 

(-8(A,)'m + 16A,(— Mfc) + 8AA,Mfe - 4(\/_ + y+)A,Mfc - —u^ 

- ^^'^ + ^(^- + ^+)«'^ + 2(rV::(r) + ry|(r))A,«fc 

- 2^{y'_{r) + V:[(r))n, - 2Ar(K:(r) + V'^{r))uk 
+ 2rV::(r)y+(r)ufe + 2rV.{r)V'_^(r)uk 

+ (rf + 2r|-)((V_ - r+)A,, - A,(r_ - Mfc), 
= 0. 

This imphes: 

[-8(A,t.,)2 - (16_ + 8A)(M,)2-(— + — Kj^'^-Mr 



+ y [4(V_ + V+) + 2r(Vl), - 6r(V+),](K)r) V-Mr 

+ '''~^''''~^'~^°" ((V+), ) - ((i - 2)(l/_)„ + 3d(l/+)„ 

r 

+ 2r((\/+),r_ + V.{y^)r) + 2r(V+),,, - 2Ar((Vl), + {y^\)\uY-^dr 
< 0. 

In the preceding calculations, we integrated by parts several times above. To justify this, 
r = must be a root of Uk of sufficiently high multiplicity to compensate for the singular 
terms. Fortunately, spherical harmonics result from eigenvalues of the spherical Laplacian. 
These take the values 



Uk = + {d - 2)k, 

and for each k, the eigenfunctions (and hence the spherical harmonics) are traces of har- 
monic polynomials of degree fc. As a result, in order to give a smooth solution as guaranteed 
above, Ufc(O) must be a zero of multiphcity k, or u^^\o) = for all m = 0, 1, ... , k. Hence, 
for k > max{0,5 — d}, the behavior of u is sufficient to make the calculations rigorous. 
See [38], Chapter 8 for a detailed description of eigenfunctions for the Laplacian on the 
sphere. 

In the commutator expression, the parameter that must dominate is a"^. Since V_, 
are smooth, exponentially decaying functions by assumption, all terms involving V-, V+ 
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and derivatives thereof are nicely bounded at and exponentially decaying. Hence, for 

4 

< r < 1, all of the functions above are easily controlled by ^ for ot sufficiently large. 
Similarly, for r > r^, sufficiently large, the exponential decay of V^_, and their 

4 

derivatives imply that any function above is dominated ^ once ot is sufficiently. In the 
intermediate region, using the smoothness of the potential functions we can find large 
enough to bound the lower order terms. In order to determine a exactly, a careful analysis 
must be done involving all of the extrema of V^_, Vj^ and their derivatives. As these 
functions are uniformly bounded, there exists such that for all a > Oq, the operator 
when conjugated by the radial Mourre operator gives a negative definite system. Hence, 
the result holds. □ 

Remark A.l. Any embedded eigenvalue can he expressed purely as a finite sum of spherical 
harmonics with radial coefficients. Combining this with the result in Section A.l\ gives a 



limited range for the calculations one must do in order to determine whether or not an 
operator has no embedded eigenvalues. 

Appendix B. Proof of the Invertibility of the Operators 



In this section we give a full proof of Proposition 6.4 This relies on our numerically com 



puted indexes, from Proposition 5.3 This proof generalizes to the 3d problem, establishing 



Proposition 5.6 



Following 13 ,20 , we first prove uniqueness and then existence. Before beginning the 
uniqueness part, we recall the following extension of the Levinson Theorem, [7 17 , from 11 

Theorem 12 (Eastham). For the equation 

+ ai(x)y("-i) + . . . + ar,{x)y{x) = 0, 

assume 

poo 

x^~^\aj{x)\dx < oo 



for some a > and j = 1, . . .n. Then there exist n solutions, yk{x), such that as x oo, 

= o(x(^-*)), k + l<i<k. 

To prove uniqueness, let u G L°°(]R) solve Cu = 0. We prove u = 0. The equation, 
—u" + V{x)u = 0, satisfies the hypotheses of the above theorem, so there exist two solutions, 
Pi{x) and P2{x), such that as x — > +oo, 

(B.l) pi(x) ~ 1, p[{x) = o(x"^), 

(B.2) P2(x)~x, p'2(x) = o(l). 
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Due to the behaviour as x — > +00, pi and p2 are hnearly independent. These can then be 
extended to all of R by the classical theory of hnear systems with smooth coefficients. 

By the same argument, there exist pi and p2 satisfying 

(B.3) ~Pi{x)r.l, p[{x)^oix-'), 

(B.4) P2ix)^x, p'^ix)^oil) 

as X ^ —00. These two are also linearly independent. We will make use of these two sets 
of functions in what follows. An important relation between them is that, by uniqueness, 
Pji^) ^ Pji-x) for J = 1,2. 

Given the values u{0) — uq and u'{0) — Uq, we know from the existence and uniqueness 
of solutions to linear systems with smooth coefficients, that there exist two unique pairs of 
constants, {01,02} and {ci,C2} such that 

u{x) = cipi(x) + C2P2{x) 
= CiPl{x) +C2p2{x). 

As x — >■ +00, 

u{x) ~ Ci + C2X. 

Since u{x) G L°°, we conclude 02 = 0. Since it is proportional to pi{x), u'{x) vanishes as 
X — >■ +00. An analogous argument ensures that u'{x) also vanishes as x — >■ — 00. Thus we 
have 

(B.5) lim u'{x) = 0. 

|a;|— >-oo 

Since u e C^(M) and its derivative vanishes, we conclude u' e L°°(R). Furthermore, we 
claim u' e L^(]R). Multiplying the equation by u and integrating by parts, 

L fL 



\u'Ydy-uu'\'l+ i Vu^dy = 0. 

I-L Jl 

Sending L — > 00 proves the claim. Additionally, this shows {Cu, u) 
Let Xa{x) be the cutoff function 



(B.6) xa{x) 



and assume A> 1. Let ua{x) — xa{x)u{x). We prove 

(B.7) lim {Cua,ua) = {Cu,u) 

A-^00 
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which will be essential to the uniqueness proof. Trivially, the potential component, J Vu\, 
will converge to / Vu"^, since V is highly localized. We now justify the convergence of the 
kinetic component. 



(B.8) 



J\u'A\dy = j XAW\'^dy + j \xA\'^\u\'^dy + 2 j XAUXA^dy 



h + h + h 



The integral Ji converges to /Im'P, the desired quantity. We must show the other two 
vanish. First, we split up I2 into 



JA 

Using our explicit characterization of the cutoff function 

^^^^^ a;(log|x|)2' 

^\ 12^ ^11 112 4(logA)^ ^ 
\Xa\ H dy < \ u\ l^ / — ^dy. 

As A —7- +00, the integral is ~ l/((log(y4~-'^))^y4), which vanishes as A ^ +00. The integral 
over {—A"^, —A) is treated similarly. 

The other integral is 

XAu'xA'^dy + 2 / XAu'xA'^dy. 

-A^ J A 

Again, using the explicit characterization of the cutoff function, 

f^' , , , ^ ,|l 21ogA ^ 

/ XAUXA'^dy <\\u\\loo\\u\\l^(a,A^) —r, — ri\2"^ 

J A J A y[^og\y\y 

— \\^\\l°°\\u\\l<^[a,a'^)- 



Since the derivative vanishes as x — > +00, (B.5), this also vanishes. The other part of Is is 
treated analogously. This proves convergence of the bilinear form. 

We now specialize to either even or odd functions. By our index computations in Propo- 
C^^^°'^ each have index 1, except for which has index 0. Without loss of 



5.3 



sition 

generality, we assume C is an index 1 operator, and proceed. Let ip be the negative eigen- 
vector with ll^/'llia = 1 for the relevant symmetry. Let Va = span j?/', m^}. We now show 
that C restricted to this subspace is negative definite. By the index computations this will 
imply ua and ip are collinear, allowing us to conclude that u = since ua has the same 
symmetry properties to u. 

Let q be any element of Va, 

q = CiUA + C2^. 
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Then, 

q) = cj {Cua, ua) + 2ciC2 {Cua, q) + cl (Cip, 

We claim this is negative, which shows C, restricted to Va, is negative definite. By the 
index computation, there is only one negative eigenvalue. Thus, dimV^ = 1, and we must 
have Ua = c{A)ip. But then 

-Xc{A) = {ua,X'4') = {uA.Cip) = {CuA,i') ■ 
Since the right hand side vanishes as A ^ oo, we conclude that c{A) = 0; hence u = 0. 
To prove the claim that the form is negative, it is equivalent to show 



(B.9) 

As A oo, 



< {Cua,ua) {C^j, 



{Cua, Ua) -5o j e~ 
{CuA,ip) -So j e~ 



\u\ 



and 



{C^,tlj) < A < 0. 



Thus, (B.9) holds for A sufficiently large and sufficiently small. Indeed, given u ^ 0, 
and ip, let 

Sn < 



1 |A|/e- 


-^y^\u\^dy 


2 Je-\y\ 


\u\\ip\dy 



Fixing this value of Sq, we can then find a value of A sufficiently large such that the 
inequality holds. 

For the operator C^°\ the proof is simpler, as we need only observe that 



ua,uaJ < 0, 

contradicting the positivity of the operator. This conludes our proof of uniqueness of the 
solutions. 



We now prove existence. Again, this follows 13,20 . We have the two fundamental sets 
of solutions {pi, P2} and {pi,p2}- Pi and pi are asymptotically constant at 00 and —00, 
respectively. Note that these two must be linearly independent, for if they were collinear, 
we would have a solution in L°°, solving Cu = 0. Hence, 

|pi(a;)| < K\x\ as X — )■ —00, 

|pi(x)| < K\x\ as X — > +00 

for some constant K. We construct a Green's function from these two to get the solution 



u{x) = pi{x) 



W{s) 



ds + pi{x) 



W{s) 



ds, 
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where W = pip[ — p[pi is the Wronskian. The integrals converge and have the appropriate 
decay due to the properties of pi and pi, and our assumption that / is highly localized. 

Finally, if / is even, consider 

= '^^-"^ L Wis) + I Wis) 

Thus = u(— x). An analogous proof holds for / odd. 



Appendix C. Numerical Methods 
The software tools we use in our computations are the Matlab and Fortran 90/95 im- 



plementations of an adaptive nonlinear collocation algorithm discussed in 26 28]. Though 
they are quite similar, we found the Fortran implementation to be faster and more robust 
for solving the 3d problems which require us to compute the ground state. For the Id 
problems, where we have an explicit formula, the Matlab algorithm, bvp4c sufficed. We 
use these tools to solve for the soliton, compute the index functions, and solve the relevant 
boundary value problems and associated inner products. 



The codes used to perform these computations are available at http: //www. math, toront 3 



edu/ simpson/f iles/ spec_prop_code . tgz^ 



C.l. Singularities at the Origin. A useful feature of this algorithm is that it can handle 
boundary value problems of the form 

y-y = -'^y + f(r,y), 
ar r 

where S is some constant coefficient matrix. The singularity naturally appears in the 
3D problems due to the Laplacian. The higher harmonics introduce a singularity which 
can be addressed by a change of variables. Let C^''^ denote one of the operators applied to 
the fc-th harmonic, 

^(fc) _ d-1 d ^ kik + d-2) 

dr"^ r dr 

Then if W solves C^'^^W = /, where / could be zero and W is non singular at the origin, 
then 

\imr~^Wir) 

r^Q 
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is a nonzero constant. This motivates the change of variable W{r) = r^W{r). In terms of 
W{r), the equation becomes 



where 



£W = _ d-l + 2k d ^ ^ 

dr^ r dr 



and W satisfies the condition 1^(0) = 0. We compute with C^''^ to get W and then muhiply 
by r''. When we compute indexes for these higher harmonic operators, the other initial 
condition becomes W{0) = 1. 

C.2. Artificial Boundary Conditions. Another subtlety of the computations is the far 
field boundary conditions. The soliton, R, vanishes as r — )■ +oo, but we only compute out 
to some finite value, rmax- For simplicity, we use the notation r and r^^x for both ID and 
3D. To accommodate this, we introduce an artificial boundary condition at rmax, and then 
check a postiori, that it is consistent. Thus, we must do some asymptotic analysis. 

We seek an asymptotic expansion for R, using (2.1) 

{-A + X- f{R))R = 0. 
As r — )■ oo, we look for an expansion of the form 



(C.l) e-^W^'Y, 



oo 

CnT 

n=0 



To extract the leading order behavior, we wish to find 7. To this end, we have 

— (e-^V^) = Ae-^V^ - 2X-ye-'^W"'-^ + 7(7 - ne'^V^-^. 



Plugging (C.l) into (2.1), we see 

r 

Hence, 7 = —{d — l)/2 and the leading order behavior is 

(C.2) r-'^e-^^ 
Therefore, as r — )• 00, 

'r^c-^'' d = l, 
i?,ie-^" d = 3. 



(C.3) Rir) 
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(a) 3d Cubic NLS (b) 3d Cubic NLS Decay Rate 

Figure 13. The numerically computed ground state for the different prob- 
lems. Computed on the indicated domain, with a tolerance of 10^^^ while 
assessing the indexes of the operators restricted to even functions for = 1 
and the zeroth harmonic for = 3. The computed solitons are monotonic 
and decay at the anticipated rate. 



From this, we develop the Robin boundary condition. 



lim ^ . , 
r-5>oo R'[r) 



d=l. 



1 

1+VAr- ' 



which we formulate as 
(C.4a) R{r 



, , —B'(r ) = for (i = 1 

max J I /-^-'I'Vinaxj " luiu, x 



and 
(C.4b) 



+ 



-B'ir 



for = 3 



assuming we have taken rmax sufficiently large. For our computations, aside from noting 
that the solver algorithm ends without errors, we have two a postiori checks available. The 
first is to verify that we have, in fact, computed the ground state. Plotting the computed R 
on both a linear and a log scale in Figure[T3]we verify that i? is a hump shaped monotonically 
decaying function. It also has the anticipated r~^e~^ decay rate. 

The second thing that can be checked is that the numerically computed R satisfies, 
asymptotically, 
-R'i 



the artificial boundary condition (C.4b). 
and observe that it vanishes as r — )■ oo. As can be seen in Figure 



To do this, 



we plot R{r) + 
the 



14 



mismatch in the artificial boundary condition is small and monotonically decaying m r 
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Figure 14. As r oo, asymptotically satisfies -AR + Ai? ~ 0. The 
figures indicate that the relevant domain, which is different for different prob- 
lems, is sufficiently large that artificial boundary conditions are good approx- 
imations. From the same computation as in Figure [13 



To compute the off axis eigenstates (p for the 3D cubic problem with eigenvalue a > 0, 
we rely on the relationship 

Asymptotically, this is the free equation 

(-A + l)Vi = -^Vi- 

Seeking a radially symmetric solution and using similar expansion techniques as in the case 
of the soliton, we will find that 

(C.5) 01 (r) ^ Cir~{d~l)/2^-rVl+ia + ^^^^{d-l)/2^-rVl=i^ 

as r — )■ oo. 
Let 

(C.6) 6 = -arctan((T), 

(C.7) p=il + aY\ 

In ID, we can construct an artificial boundary condition that, as r — )■ rmax, 

(C.8) 0; + pcos(0)0i+psin(^)02 = 0, 

(C.9) 02 + pcos(^)02 -psin(^)0i = 0. 
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(a) Unstable Mode of 3d Cubic NLS (b) Unstable Mode 3d Cubic NLS Decay 

Rate 

Figure 15. The numerically computed ground state and the off axis un- 
stable mode of 3d cubic NLS. Computed on the indicated domain, with a 
tolerance of 10~^^. 



Analogously, in 3D, as r — )■ rmax, 
(C.IO) (t)\ + p cos 

(C.ll) 02 + pCOS 



1 + -</>i + psin(^)(/)2 = 0, 

r 

2 + -02 - psin(6')0i = 0. 



As in the case of the soliton, we can verify that the solutions have the appropriate shape, 
decay as expected, and satisfy the artificial boundary conditions. The shape and decay are 



plotted in Figure [T5j The functions rapidly reach machine precision. If we zoom in on the 
unstable modes, as in Figure [16| we can see the periodic structure. However, as is suggested 
by these figures, once 0j's are sufficiently small, < 0(10^^^), this fine structure degrades. 
Fortunately, this numerical error is sufficiently small as to not impact our computations. 



The artificial boundary condition plot appears in Figure 17 



Another place where we use artificial boundary conditions is in solving the various bound- 
ary value problems for t/j'^^ or Z^^\ Tthe rapid decay of the soliton leads to the functions 
satisfying the free equation, for (i > 1, 



- Ag + k{k + d-2) /r''-'+^q = 0, 

where q is one of these functions. In this region, the function must asymptotically be like 
(5.7b). Of course, we work with the variables f/-'^^ and zf \ Since these vanish as r — )■ oo, 



they asymptotically behave as 



q (X r 



2-d-2k 
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Figure 16. Zooming in on the components of the numerically computed 
unstable mode, 0, we see the periodic structure of the solution. 



Thus, we have the artificial boundary conditions, valid for all harmonics and d > 2, 
(C.12a) f/r(rmax) + 



(C.12b) 



^ v n 



+ 



d-2 + 2kdr ^ 
r d 

d-2 + 2kd^"^ 



0. 



Figure [17] shows that these artificial boundary conditions are asymptotically satisfied. 
For our d = \ computations, we have the artificial boundary conditions 



(C.13) 



d 

dx 



d 



0' 



0. 



We can similarly check that the inner products are asymptotically constant and that our 
computed functions asymptotically satisfy the artificial boundary conditions. See Figures 
[Island [021 
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Figure 17. As r — )• oo, our functions satisfy the asymptotically satisfy 
the free, linear equations for the different parts of the 3d cubic problem. 
The figures indicate that the relevant domain, which is different for different 
problems, is sufficiently large that artificial boundary conditions are good 
approximations. 
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(a) 



C^' with a = 2.1. 
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(b) /:^°^with (7 = 2.1. 
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(c) with a = 2.5. 
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(d) £^°Vith a = 2.5. 
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(e) C^^^ with cr = 3.0. 
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(f) = 3.0. 



Figure 18. As r — )■ oo, our functions asymptotically satisfy the free equa- 
tions for the various Id supercritical problems. 
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(a) Cj.'^ in the critical case with the (b) /I^-' in the critical case 
natural orthogonality conditions. 
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(c) in the critical case with (d) in the critical case with the 

the alternative orthogonality con- FMR orthogonality conditions, 
dition. 

Figure 19. As r — )• oo, our functions asymptotically satisfy the free equa- 
tions for the Id critical problem with various orthogonality conditions. 



The last place we make use of artificial boundary conditions is in solving = —x^R 
for the critical problem in Section 6.2.3[ Using the same procedure as above, one will find 
that (3 oc x^e~^ form which the artificial boundary condition 



(C.14) 

can be constructed. 



/3'(r) + ( 1 - ^ ) m 



C.3. Computation of the Indexes. In computing the indexes of operators C^^^ (from 
which we recover the indexes of C^±''), we simultaneously solve the mixed boundary value 
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problem/initial value problems 



(C.15) -AR + XR- g{\R\')R = 0, R'{0) = 0, (C.4a) or (C.4b), 



(C.16) rPf/^'^) = 0, U^''\0) = 1, ^U^''\0) = 0, 

dr 

(C.17) = 0, Z^''\0) = 1, ^Z('')(0) = 0. 

dr 

For the 3d cubic equation, this is the complete set of equations; A = 1 and f{s) = s. 
The analogous computations are made in ci = 1. In addition to verifying that the soliton 
was adequately computed, we can check, a postiori, that the index functions, f/^^^ and 
Z^^^ asymptotically satisfy the free equation. This was the shown in the index figures of 



Sections |5.2| and 6.1 where we checked the constants. 

C.4. Computation of the Inner Products. As in the computation of the indexes, we 

/j") and 



similarly solve mixed boundary value/initial value problem for R, the ?7-"^ and Z^°'\ and 



the inner products 

In computing the inner products, we introduce the dependent variables K^°'\r) and 
7j-"'*(r), where 



(C.18) ^nfir) = C^M:Hfr'-\ ^f{0) = 0, 

d 



dr 



(C.19) (r) = Ci"^ Z^f Z;,:^\''\ 7]"^(0) = 

for ii and £2 the appropriate indexes. Clearly, 

lim Kf{r) = K^''\ 

lim 7,-°^(r) = 



We approximate the inner products by computing to rmax- As demonstrated in Figures 20 



21 , 22 these converge rapidly and are essentially constant in the region of the free equation. 



This is entirely consistent with the exponential decay of the soliton and functions related 



to it, such as its derivative. 
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(a) Inner Products for C^^^ functions for 3D (b) Inner Products for C^^'' functions for 3D 
Cubic Cubic 
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Figure 20. Inner products for the 3d cubic equation. 
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Figure 21. Inner products for the Id supercritical problems 
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(c) £^'' in the critical case with (d) C'"_^' in the critical case with the 
the alternative orthogonality con- FMR orthogonality conditions, 
dition. 
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Figure 22. Inner products for the Id critical problem with various orthog- 
onality conditions. 
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